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<D Abstract 

g 

We derive an effective Bose-Hubbard model that predicts a phase transition from Bose-Einstein 
condensate to Mott insulator in two different systems subject to applied periodic potentials: mi- 
crocavity exciton-polaritons and indirect excitons. Starting from a microscopic Hamiltonian of 
electrons and holes, we derive an effective Bose-Hubbard model for both systems and evaluate 
the on-site Coulomb interaction U and hopping transition amplitudes t. Experimental parameters 
J> required for observing a phase transition between a Bose-Einstein condensate and a Mott insulator 

o 

^sO are discussed. Our results suggest that strong periodic potentials and polaritons with a very large 

excitonic component are required for observing the phase transition. The form of the indirect 

cn 

exciton interaction is derived including direct and exchange components of the Coulomb interac- 

o 

t-H tion. For indirect excitons, the system crosses over from a Bose-Hubbard model into a double 

• — h layer Fermi-Hubbard model as a function of increasing bilayer separation. The Fermi-Hubbard 

X 

model parameters are calculated, and the criteria for the location of this crossover are derived. We 
conjecture that a crossover between a Bose Mott insulator to a Fermi Mott insulator should occur 
with increasing bilayer separation. 

PACS numbers: 71.36.+C, 71.35.-y, 03.67.Ac 
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I. INTRODUCTION 



The observation of the Bose-Einstein condensation (BEC) of exciton-polaritons has gen- 
erated a large amount of interest in recent years [IH3] • The focus has now turned to exam- 
ining various properties of the condensate, such as thermal equilibration jl], superfluidity 
[5], vortex formation j6j [7], and elementary excitations [8]. If the trend followed by atom 
optics physics community holds for the exciton-polariton community, one important branch 
of study of exciton-polariton BECs will be the application of periodic potentials on the 
BEC system. Optical lattices have attracted much attention, spurred on by the experiments 
demonstrating the phase transition between a BEC and a Mott insulating state in a Bose- 
Hubbard model [9]. The application of the periodic potential simultaneously increases the 
particle-particle interaction, as well as decreasing the kinetic energy. This allows the ratio 
of the Hubbard on-site interaction to the hopping amplitude U/t to be varied at will. The 
experiment has been of particular interest in the quantum information community, since 
the experiment realizes a nearly ideal quantum simulator [TUl ITT] . A quantum simulator 
is a device that directly recreates a quantum many-body problem in the laboratory. By 
experimentally modifying physical parameters, such as the periodic potential amplitude, 
temperature, and density, one may explore the phase diagram of the system. 

The formation of a Bose-Hubbard model using polaritonic systems was first proposed in 
refs. [121 H3]. In this paper we develop the theory for exciton-polaritons subject to a periodic 
potential (see Fig. [T^). In contrast to the works of refs. [T31 H31 113 HHHZJJ, where the 
polariton interaction originates from an effective nonlinearity due to a coupling to atomic 
sites, our interaction originates from the excitonic components of the polaritons, which 
ultimately originates from a Coulomb interaction. Starting from a microscopic Hamiltonian 
for electrons and holes and their Coulomb interaction, we derive the origin of the Bose- 
Hubbard model that is assumed in Ref. [IB], allowing an accurate determination of the 
Bose-Hubbard parameters U and t. From an experimental point of view, steps towards 
a similar experimental configuration as the optical lattice have been realized already by 
modifying the semiconductor microcavity system. In ref. [H], it was shown that a band 
structure was successfully formed using a metal deposition technique. The periodic metal 
structure on the surface changes the boundary conditions of the photon field, thus creating a 
static periodic potential for the polaritons [T3] . Other methods for trapping polaritons have 
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been proposed by etching the microcavity [IB] . Such etching techniques are anticipated to 
produce stronger trapping potentials and access a more strongly correlated regime [17] . 

The formal similarity in the treatment of indirect excitons allows us to write general 
formulas that capture both the polariton and indirect exciton interaction (see Fig. [T|b) . 
Polaritons are described in the d = limit of the formulas, where d is the bilayer separation 
of the indirect exciton system. Indirect excitons have a non-zero d, but a vanishing photon 
component. Although only exciton-polaritons have currently been observed to undergo Bose- 
Einstein condensation so far [221123], there is a large amount of interest in BECs of indirect 
excitons, as well as works of indirect excitons in periodic lattices [24], motivating us to write 
the generalized formulas for both cases. We place particular interest on what parameters 
are required for observing a Bose-Hubbard Mott transition. For indirect excitons, as the 
bilayer separation d is increased, the bosonic nature of the excitons gradually diminishes 
due to the reduced electron-hole interaction. The system is more appropriately described as 
a Fermi-Hubbard model in this limit. We discuss the criterion for this crossover to occur for 
our model. We also conjecture that a crossover between a Bose Mott insulator to a Fermi 
Mott insulator should take place with increasing bilayer separation, in analogy to the more 
commonly known BCS-BEC crossover [39] 140]. 

SI units are used throughout this paper. 

II. BOSE-HUBBARD MODEL 

We assume that a periodic potential of the form 

W ph (r) = Wt [cos(k x) + cos(k y)} 



is applied on the photonic (ph) and excitonic (exc) components of the exciton-polaritons 
respectively, where ko = 2tt/\ and A is the wavelength of the periodic potential created. 
As mentioned in the introduction, a variety of experimental methods exist to create such 
a potential on the photon field [I4T4IB] 125] . For the excitonic part, metal gates may be 
applied to the surface, trapping the excitons under the gates [2BJ. The potential W eKC { r ) ls 
an effective potential for the center of mass motion of the exciton, obtained after integrating 
over the relative motion of the exciton. The potential can in principle be either type I (where 




(1) 
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the electron and holes share the same potential minimum locations) or type II (electrons and 
holes have minima on alternate sublattices). For example, deformation potentials induced 
by surface acoustic waves [2Z] and dipolar traps J2S] are type I potentials. Meanwhile, the 
piezoelectric trapping technique [28] is an example of a type II trapping potential. Type 
II potentials need rather strong trapping potentials for each individual component of the 
exciton (electron and hole) since their effective amplitude is suppressed by a factor of (/coos) 2 
(see. eq. (30) in Ref. [28]). However, since for strong potentials type II potentials tend to 
ionize the excitons, we believe that a type I potential is more promising in order to avoid 
these undesired effects. 

The total Hamiltonian of the system is then 



H = H cxc + H c , 



H cxc = / d 2 rb\ 



~\~ H ph + -ffexc-ph ~t~ H, 

h 2 



6(r) 



H, 



P h 



2M 

d 2 Qd 2 Q'd 2 qtf Q _ q tf QI+q UUQ, Q', q)b Q >b Q , 
h 2 
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2m 
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H. 



cxc-ph 
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where 



b(r) 



air) 
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2^ 
1 

2^ 



d 2 qe lk - r b 



d 2 qe iqr a q 



(2) 
(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
(9) 



are the annihilation operators for the quantum well excitons and microcavity photons respec- 
tively, M = m e + rrih is the exciton mass, g is the exciton-photon coupling, and U cxc (r,r') 
is the effective interaction between two excitons. The photon acquires an effective mass 
m ph through the dispersion in a two-dimensional microcavity, where the photon energy is 
E ph = m ph c 2 , and c is the speed of light in GaAs. We do not consider the spin of the excitons 
explicitly because we assume that the polaritons are injected with a linear polarization such 
that only one spin species is present. The form of the exciton interaction is discussed in 



detail in section HI A H mt is the non-linear interaction due to the exciton-photon coupling 



1321 and is discussed in section III B 
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Substituting eqs. ([§]) and ^ into (J2|, and denning the upper (cr =|) and lower (a =!) 
polariton operators 

^ = u a b q + v a a q , (10) 
where \u a \ 2 + |f°"| 2 = 1. We obtain the polariton Hamiltonian 



H = H kin + H pot + H pol 



H, 



pol 



2 ,t, 

\u\ 2 W^ + \v\ 2 W» h ) f J2 r t t 
d 2 Qd 2 Q'd 2 qU pol (Q, Q', q)PQ„ q PQ /+q PQ>PQ 



u POl (Q,Q',q) 



\u\ 



-u eKC (Q,Q',q) 



\u\ 2 (uv* 



U V 



uUQ,Q',Q-v) 



(ii) 

(12) 
(13) 
(14) 
(15) 



2 - cxc \ -v : -v ; 1/ 2 

where we have only included terms where lower polariton operators appear, and compacted 
the notation such that p q = p q . Physically, disregarding the upper polariton operators cor- 
responds to a low-temperature regime where there is negligible upper polariton population, 
which is routinely achieved experimentally. The Hopfield coefficients u = and v = ir are 
taken around q = 0, again assuming that only the low energy states are excited. The lower 
polariton dispersion is obtained by expanding around q = giving 

h 2 q 2 



2m 



where the lower polariton mass is given by 



m 



P ° ' 



\u\ 



pol 



+ 



rn 



P h 



Reverting to real space makes it clear that we have polaritons in a periodic potential 

h 2 



H 



2m 



-V 2 + \v\ 2 W pb (r) + \u\ 2 WUr) 



pol 



p(r) + #, 



pol * 



(16) 



(17) 



(18) 



For sufficiently low temperatures, we may retain the lowest energy band of the Hamilto- 



nian (18) to a good approximation. A necessary temperature criterion is that the thermal 



energy /cgT is less than the band gap A. In one dimension, any non-zero potential Wq will 
open a bandgap. In two dimensions however, for small potentials the lowest energy band 
overlaps in energy with the second lowest energy band. To ensure that these bands are 



separated, a potential of approximately Wo 

W = \v\ 2 W^ h + |4 2 Wo e 



— is needed 1291. where 



2m pol 
2TJ/Ph I I „ . I 2 TT ■ ' 



(19) 



is the total potential amplitude due to exciton and photon parts. Under these circumstances, 
we may make a Wannier transformation and retain only states in the lowest energy band. 
This yields 

H = t(n, n')p\j) n i + - U (ni, n 2 , n 3 , n^jp^p^p^p^ (20) 



ni,ri2,n 3 ,rL4, 



where 



t(n,n') = J d 2 rw*(r — n\) 



h 2 



2m 



-V 2 + \v\ 2 W ph (r) + \u\ 2 WUr) 



pol 



w(r-n'\), (21) 



{/(ni, 712,713,714) = J d 2 rd 2 r'w*(r — ni\)w*(r' — n 2 \)U pol {r, r') 

xw(r' — n 3 X)w(r — n 4 A), (22) 



and 



p n = I d 2 rw(r — n\)p(r). (23) 
w(r — nX) is the Wannier function centered around the lattice point n = (n x ,n y ). The 



Hamiltonian (20) is a Bose-Hubbard Hamiltonian. We shall be only concerned with the 



nearest neighbor tunneling matrix elements t and the on-site Coulomb interaction U in this 
paper: 



tS ( <».» + <i,o)) =/«**•<*) 



h 2 



2m 



- V 2 + W cos(k x) 



pol 



w{x-X), (24) 



and 



U = U (n, n, n, n) 



A|2 



d 2 rd 2 r'\w{r)\ 2 U pol {r,r')\w{r')\ 



(25) 



In writing (20) we assume that the polariton lifetime r is sufficiently long such that the 



extended "superfluid" and Mott states can occur. For example, for an extended "superfluid" 
state we require that there is enough time for the polariton to hop several times before 
decaying 



T 



> 



h 
1 



Similarly for the Mott insulating state, the Coulomb energy should obey 

h 



T > -. 

~ u 



(26) 



(27) 



In addition to the bandgap criterion ksT < A, it is also necessary to have /c#T < U, t, such 



that only the low energy physics of the Bose-Hubbard model is probed. To evaluate (25) we 



require a form for the polariton-polariton interaction which we evaluate in the next section. 



III. EFFECTIVE POLARITON INTERACTION 



A. Exciton-exciton interaction contribution 

To obtain the effective interaction between excitons, we use the methods of de-Leon and 
Laikhtman [30] . There are many approaches to find the effective interaction between excitons 
in the literature, such as usage of the Usui transformation [31], variational wavefunction 
methods [33], and operator methods [34]. We find that the wavefunction methods of Ref. 
[50] are most transparent and systematically give the exciton interactions to order a B /A, 
where a B = Aireh 2 /2e 2 fi is the 2D Bohr radius, A is the trapping area of the excitons, and 
fj, is the reduced mass /i = m e m,h/{m e + rah). We note that a similar method was used by 
Ciuti et al. [35] to find the same result for the Coulomb interactions, but it is unclear how 
to treat corrections to the kinetic energy operator ( "kinematic corrections" ) based solely on 
their method. Ref. [30] makes it clear that such terms cancel in the end and do not give rise 
to a physical interaction. Here, we generalize the results to indirect excitons in a periodic 
potential. 

Following Ref. [30], the effective Hamiltonian for the excitons can be decomposed into 
the following terms (omitting function labels for brevity): 

Uexc =U dit + U* ch + Ul ch + Ut h - 1 -{AU + UoA) + (28) 

where the terms are the direct exciton scattering, the exciton exchange scattering, the 
electron exchange scattering, the hole exchange scattering, the correction due to non- 
orthonormality of the wavef unctions, and the contribution due to excited states of the 
excitons. Since only exciton-exciton interactions to order a B /A are kept, only two-body 
exciton interactions need to be considered. Explicit expressions for the above terms are as 
follows. The direct term is 

U diI (Q,Q',q) = J rfV e rfV /l rfV e /rfV /l /^Q(r e ,r /l )^Q,(r e /,r /l /)^^ Q+(? (r e ,r /l )^Q/_q(r e /,r /l /), 

(29) 
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where Hamiltonian for the two exciton system is 



"« = ~ - - 4v?„ - V(W - n\) - v(W, - r, A) 

Hl = V (\ re _ re ,\) + V (\r h - r h ,\) - V(^/\r e - r h ,\ 2 + d 2 ) - V(y/\r h - r e ,\ 2 + d 2 ) (30) 
with V(r) = e 2 /47rer (e ~ 13eo is the permittivity in GaAs). The exciton exchange term is 

UZJQ, Q', q) = U diz (Q, Q', Q'-Q-g). (31) 



The electron (hole) exchange terms are obtained by multiplying (29) by —1 and exchanging 
r e -H- r' e (rh -H- r' h ) in the final states. The correction due to non-orthonormality is 

A(Q, Q', q) = -I e A (Q, Q', q) - /J(Q, Q\ q) (32) 

where 

Il(Q,Q',q) = J d 2 r e d 2 r h d 2 r e ,dh h ,y* Q (r e ,,r h )^,(r e ,r h ,)^ Q+q (r e ,^ 

I h A {Q,Q'iq) = J d 2 r e d 2 r h d 2 r e ,d 2 r h ,y* Q (r e ,r h ,)^* Q ,(r e , ,r fc )* Q +q {r e ,r h )^Q^ q {r e ,,r h/ ). 

(33) 

The final term l{( exex ) originating from excited states of the excitons is neglected in our 
analysis since its magnitude is small, using similar arguments as given in Ref. [30J. 

To evaluate the expressions above, we use an approximate form for the indirect exciton 
ground state wavefunction, as obtained in ref. [37] : 

V Q (r e ,r h ) = -^=exp[iQ ■ ((3 e r e + f3 h r h )]G(p, Z)f^(z e )f^ h \z h ), (34) 
V A 



where p = ^J(x e - x h ) 2 + (y e - y h ) 2 , Z = z e - z h , and /3 e>h = m e , h /(m e + m h ). In eq. (J34J), 
the exciton is considered to be trapped in a large area A, such that the center of mass 
wavefunction is of the form of a plane wave. The relative wavefunction of the exciton is 

G(p,Z) = ^exp \-(X(Z)/2)(y/(p/a B y + (Z/a B ) 2 - Z/a B )] , (35) 
a B L J 

where X(Z) = 2/(1 + y/2Z/a B ) and N G is a normalization factor such that 
J \^Q(r e , rh)\ 2 d 3 r e d 3 rh = 1. P e \z) (P h '(z)) is the electron (hole) wavefunction satisfy- 
ing the Schrodinger equation in the ^-direction. For indirect excitons with the electrons and 
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holes confined to their respective quantum wells, we may approximate the wavefunctions to 
be f^ e \z e ) = 5(z e — d/2), f^ h \zh) = 5(z h + d/2), and Z = d. The exciton wavefunction for 
polaritons may be recovered by setting d = 0. 



The evaluations of the direct and exchange terms are deferred to Appendix A 1 We find 
the direct term to be 



U diI (Q,Q',q) 



-2Eu 



2M 



{Q 2 + Q 



I2\ 



6{q) + jh aB 



I dit (q,d). (36) 



where E\ s is the binding energy of a Is exciton. The function I dil (q, d) is plotted for various 
d in Fig. [2j The exciton exchange term is 

h 2 



U*JQ,Q',q) 



—2E\ S 



2M 



(Q 2 + Q 



I2\ 



8{Q -Q' + q) 



i 2 / n \ 2 

+ J^- e a B[-) UW{AQ) 2 + q 2 - 2AQqcos9,d), (37) 

where AQ = \Q' — Q\ and 9 is the angle between Q' — Q and q. The electron and hole 
exchange terms are 



U^ h (Q,Q' iq ) 



-2E ls + ^(Q 2 + Q' 2 ) + ^((Q + q) 2 + (Q' - q) 2 ) 



1 e 2 



A Aue 



a B 



(AQ,q,9,P e ,d) 



utAQ,Q',q) 



~ 2Els + §i {q2 + Q ' 2) + \m m + 9)2 + {Q ' ~ q)2) 



AAixe 



a B 



-^exch 



I e A (Q,Q',q) 

(38) 

Ia(Q: Q'> q) 

(39) 



where in the square brackets we evaluated half of the operator on the initial states and half 
on the final states. Doing this we see that these terms exactly cancel with the corrections 



due to non-orthonomality (i.e. the fifth term in (28)). Numerical evaluations of the exchange 
integral J cxch are shown in Fig. |3j 



Substituting (36 39) into (28) we obtain the final effective Hamiltonian for the two-exciton 
system. Subtracting the kinetic energy and binding energy terms, we obtain an expression 
for exciton-exciton interaction 



uUQ,Q',q) 



A Aire 



a B 



71 



J dir (g, d) + I dix ( y /(AQy + q*-2AQqcos6, d) 



-J cxch (Ag, q, 9, (3 e , d) - J cxch (AQ, q, 9, (3 h , d) 



(40) 
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The dependence on d for small q may be evaluated exactly for the direct term by expand- 



ing the term in the square brackets in (Al). We obtain 

de 2 



U Ait (Q,Q',q = 0) 



eA 



(41) 



This is the zero momentum limit of the Fourier transform of the Coulomb interaction for 
oriented dipoles U(q) = e 2 (l — exp(—dq))/eq. The d dependence of the exchange term must 
be evaluated numerically, and our results are shown in Fig. |4| We find an approximately 
linear dependence of the exchange term with d, which changes sign at d/as ~ 0.66. In Fig. 
[4] we also plot the combined contributions of the direct and exchange integrals. We see that 
the total interaction remains repulsive for all d, despite the exchange term changing sign. 

B. Saturation contribution 



The saturation contribution to the polariton interaction comes from the coupling of the 
electron and holes to the electromagnetic field 



#em = HjS I d 2 kd 2 k' 



e k h k ,ak + ki + a k+k ,hk'Gk 



(42) 



The last two terms in (]2]) may be found by considering the matrix element between the 
states 

\Q,Q') = ^ J d 2 r e d 2 r h d 2 r e ,d 2 r hl ^ QQ ,(r e ,r h ,r el ,r h ,)e\r e )h\r h )e\r e ,)h j (r h/ )\0) (43) 



and 



Q", q) = a\ / d 2 r e d 2 r h V Q „(r e ,r h )e\r e )h\r h )\0) 



where 



$QQ>(r e ,r h ,r e ,,r h , 



(44) 



^ Q (r e , r h )m Q ,(r e ,, r h ,) + ^ Q (rv, r h ,)V Q ,(r e , r h ) (45) 



- m Q (r e ,r h ,)m QI (r el ,r h ) - ^ Q (r e/ ,r h )^ Q/ (r e ,r h/ 



(46) 



Starting from the two exciton wavefunction, if EM can either destroy one of the excitons and 
produce a photon, or take an electron and hole from each of the excitons and produce a 
photon. These two processes give rise to the matrix element 

(Q + Q'- q, q\H EM \Q, Qf) = hGVAG(0, 0)(5(Q' ~q) + $(Q ~ q)) 

+U mt (Q,Q',q) + U mt (Q',Q,q), (47) 
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where 

MQ, Q'. ^) = / dWr e d 2 r^ Q (r e , x)* Q ;(aj, ^q+q^K, r h )e**- m , (48) 
and we have set Q" = Q + Q' — q due to momentum conservation. The first two terms in 



(47) correspond to the destruction of an exciton to create a photon, with another exciton 
acting as a bystander. This corresponds to the second last term in (|2]). The last two terms 
correspond to an electron and hole being taken from each exciton, resulting in a new exciton 
being formed from the remaining electron and hole. This process clearly requires two excitons 
in the initial state, giving the non-linear last term in (|2]). We only consider the case of zero 



bilayer separation (d = 0) here, since the Hamiltonian (42) requires that the electron and 
hole recombine into a photon at the same spatial position. In our approximation (34) where 
the electron and hole wavefunctions are perfectly confined to their respective quantum wells, 
for non-zero bilayer separation there is zero overlap of the electron and hole wavefunction, 



which gives a zero matrix element for (47). 
We thus find 

g = gVAG(0,0) (49) 

and 

EUQ, Q', i) = -hG^=iUQ, Q', <?)• (so) 

V A. 



where J sat (Q, Q',q) is evaluated in the appendix. A numerical evaluation of the integral as 
a function of the photon momentum is shown in Fig. [3] 



IV. MOTT-HUBBARD TRANSITION 



Returning to the effective Bose- Hubbard Hamiltonian of (20), we may now estimate the 



size of the tunneling and Coulomb energies from eqs. (24) and (25). Starting from (25) we 
make a change of variables -R C m = R + R' and p = (R — R')/2, giving 



U 



d 2 R CM d 2 p\w 1 {R CM + p/2)\ 2 \ Wl {R CM - p/2)\ 2 U pol (p) 



(51) 



From Figs. [2] and [3] we see that interaction is large up to a momentum of the order of 
~ 1/ag. Thus the largest contributions to the above integral occur when the variable p is 
of the order of ~ a#, which is much smaller than the length scale of the Wannier functions 
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: M I I 



where 



A. We may thus write 

= U conl + £/ sat (52) 

U ma = |w| 4 f/ cxc (0,0,0)A J d 2 R\w 1 {R)\ A 

= -\u\ 2 {uv* +u*v)U^ t (0,0,0)A J d 2 R\w~ / (R)\ 4 . (53) 

Fig. [5] shows the two contributions to U as well as the hopping amplitude t. The results 
are normalized to the characteristic energies 

(54) 



j rcoul 
^0 


2e 2 a B \u 4 
7r 3 eA 2 


j rsat 
^0 


= 2hgJ^\u 


*0 


h 2 


8m pol A 2 ' 



§ (55) 



(56) 



where we have used (49) to convert Q into g. Here, 2hg is the Rabi splitting of the polaritons. 
For polaritons, the bilayer separation is d = 0, while for indirect excitons the exciton com- 
ponent is \u\ 2 = 1. We see that increasing the potential strength Wq decreases the hopping 
t while increasing U, as expected. Increasing d enhances ?7 coul , as the dipole moment of the 
excitons are enhanced with an increasing d. 

We now derive a criterion for a quantum phase transition from a BEC state into a Mott 
insulator state. In two dimensions, the phase transition is expected to occur at approximately 
[38] 

U/t w 23. (57) 

By turning up the potential Wo it is clear that at some point U/t will reach this critical 
amplitude. The other variable that may be changed to reach the phase transition is the 
detuning of the polaritons which changes the polariton mass. For a potential of size Wq ~ 



h 2 k? 



2m 



B -, we may derive a criterion for the polariton mass necessary to reach the phase transition 



pol 



using (57) and the ratio of the dimensionless parameters in Fig. |5j For GaAs, this gives 

Tlhh 2 



m pol w HT 3 : 7 ,. w 1(T 2 M, (58) 
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where in the second relation we assumed that \u\ 2 is of the order of unity, and a# ~ 10 nm. 
This corresponds to extremely far blue-detuned polaritons, i.e. very exciton-like polaritons. 
The lack of the dependence on the wavelength A is due to the cancellation of the dependence 
of the Coulomb interaction energy Uq and the kinetic energy to- The wavelength A is still 
important however, as it sets the energy scale of the whole Bose-Hubbard Hamiltonian. The 
energy scale should be set such that the parameters U and t are larger than the temperature 
of the experiment, such that only the lowest energy band is occupied. Furthermore, semi- 
conductor systems possess an inevitable disorder potential due to reasons such as crystal 
imperfection and damage during fabrication, thus A should be set small enough such that 

the hopping energy t overcomes this disorder potential strength. For larger potentials than 

h 2 k 2 

Wo ~ 2m ° i , a lighter polariton mass is allowable for reaching the phase transition. Thus in 
practice a combination of blue-detuned polaritons and large potentials is probably the most 
favorable experimental configuration. 

For example, using typical experimental parameters for polaritons (d = 0) in GaAs using 



the criterion (58), corresponding to u ~ 0.999, with A = 0.5 jum, a# = 10 nm, 2hg = 15 
meV, and an applied potential of Wq = 6 meV, we obtain U = 0.24 meV and t = 9 /xeV. This 
corresponds to temperatures in the vicinity of T « 0.1 K, which are reachable using today's 
refrigeration methods. In order that the system is stable in the Mott insulating state, the 



lifetime of the polaritons should be longer than the timescale set by (27). Assuming that the 
lifetime of the very exciton-like polaritons can be approximated by typical exciton lifetimes 
t « Ins [44], this corresponds to an energy scale U > h/r ~ 1/xeV. We thus see that for the 
above parameters the lifetime requirement is satisfied. 

The Coulomb interaction is increased for indirect excitons as shown in Fig. [5| However, 
the increase is fairly modest for bilayer separations of the order of the Bohr radius. Thus 
considering that indirect excitons have not been observed to undergo BEC yet, the moderate 
advantage of increased interactions (at the sacrifice of a lighter polariton mass) is outweighed 
by the difficulty of cooling the system into the ground state. 

The state of the polaritons may be measured using standard photoluminescence measure- 
ments that measure the coherence across the condensate [36J . In analogy to the experiments 
of Greiner et al. [9], the transition to the Mott insulator state should lead to a lack of spatial 
coherence across the sample, resulting in the destruction of the far-field interference pattern. 
However, the disappearance of interference fringes does not unambiguously demonstrate the 
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presence of a Mott insulator, since an uncondensed state will also have the same interference 
characteristics. Therefore, a second order coherence Hanbury-Brown-Twiss measurement is 
also necessary to determine the correlations between the photons emanating from the sam- 
ple. At unit filling in a Mott insulator state, the conditional probability of detection of a 
photon originating from a particular site following detection from the same site vanishes. 
A similar diminished probability is present at higher filling factors. This is the identical 
technique used to observe anti-bunching behavior in single photon generation. 



V. FERMI-BOSE CROSSOVER BOUNDARY 



We now turn to the effect on the particle statistics of indirect excitons as the bilayer 
separation is increased. For d = 0, excitons are well approximated by bosons for sufficiently 
low density, which motivated us to describe the system as a Bose-Hubbard model in section 
[TTj In the limit d = oo, excitons cannot be described by bosons, and are more properly 
described as a double Fermi-Hubbard model, with electron-hole interactions between them. 
Since the two descriptions are rather different with differing phase transition criteria, it is 
of interest to know at what d this crossover occurs. The bosonic nature of (or the lack of) 
the excitons may be seen by examining the commutation relation 

[b, 6 f ] = 1 - K (59) 

where 

b ] = J d 2 r e d 2 r h ct) d {r e - r h )w(P e r e + (3 h r h )e\r e )h) {r h ) (60) 

and 

K = J d 2 r e d 2 r h d 2 r(p d (r e -r h )w((3 e r e + (3 h r h ) <j> d (r e - r)w((3 e r e + (3 h r)h ] (r h )h(r) 

+M r h - r)w{/3 e r + f3 h r h )e\r e )e(r)j . (61) 

We have omitted the label n and written W n — > to simplify the notation. Only the 
commutation relation in the same potential minima of the periodic potential is considered 
here (i.e. the same Wannier function), as deviations from bosonic behavior should be most 
apparent in this case. The operator K may be interpreted as the correction operator to 



the commutation relation (59), as it contains the non-bosonic component of the exciton 
operators. 
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Due to the presence of the K operator, n-particle states defined using the operators 



defined in (60) do not have the simple l/yn\ normalization of bosonic states. We must 



instead define such states according to 

In) = - 7 =^=(&tr|0), (62) 

where N{n) is present for proper normalization. A derivation of this normalization factor is 
given in Appendix |Bj up to powers linear in the correction operator K . We obtain 

N(n) = 1 - in(n - l)I BF (W Q ,d), (63) 



where I BF (Wo,d) is the integral expression given by (B8) and has an order of magnitude 



~ a 



Bl 



/A 2 . All terms neglected in (63) have higher powers of a 2 B /X 2 , which have a small 



contribution for the typical periodic potential dimensions that are possible using current 
fabrication methods (A ^> ob). 

Defined in the way ( [62] ), the states \n) provide an orthonormal basis set. Deviations from 
bosonic behavior occur due to the operator not providing the correct mapping between 



these states (n\-j=\n — 1) ^ 1. Using the definition (62), we have 



/ , bi , ,x / N(n) 

(n\-=\n-l) = \ AT , v \, . 64 
V N n-1 ( J 



Deviations from unity of the RHS represents non-bosonic behavior. Substituting (63), this 
factor is to lowest order in I BF (W , d) 



A( "' - ' 1 (n-l)I BF (W ,d). (65) 



N(n - 1) 2 

Since the Bose to Fermi transition is a smooth crossover strictly speaking it is arbitrary 
where to mark the boundary. However, a reasonable criterion for the location of the crossover 
from bosonic to fermionic behavior may be defined as when the second term in the above 
expression becomes of the order of unity: 

^(n-l)I BF (W ,d) = 1. (66) 



The solution to the above criterion is plotted in Fig. [6] in the space (d, Wo). We see that 
with decreasing n, d and Wo, the excitons become more boson-like. This dependence on 
n and d is a restatement of the well-known result that the excitons become non-bosonic 
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when their wavefunctions start overlapping, i.e. na 2 B /A ~ 1. The dependence on Wo may be 
understood by considering the spread of the Wannier functions with W . As W is increased, 
the Wannier functions become more localized, effectively reducing the area that the excitons 
are confined in. This enhances the overlap between the excitons, thus pushing the boundary 



towards the fermion side of the crossover. For n < 1, there is no solution to (66), meaning 
that the boundary for bosonic behavior extends all the way to infinity in Wo and d. The fact 
that solutions for bosonic behavior exist with n > 1 means that in reality one cannot treat 
the excitons completely as hard- or soft-core bosons, and their true nature lies in between 
these two limits. 



VI. FERMIONIC DESCRIPTION OF THE BILAYER SYSTEM 



For parameter regions where the bosonic approximation is invalid, we must write the 
Hamiltonian in its full form involving both electron e a (r) and hole h a {r) operators: 



h 2 



2m, 



-V 2 + W e {r) 



e -( r ) + Yl J d2rh l( r ) 



h 2 



2m h 



-V 2 - WJr) 



hJr) 




d 2 rd 2 r' 



el(r)et(r')V(\r-r'\)e a (r')eAr) 



+ hl(r)hl(r')V(\r - r'\)K(r')h a ,(r) - 2el(r)hUr')V(^\r - rf + d^)K(r')e al (r) 



(67) 



where 



W e (r) = W a \cos(kox) + cos(koy)} , (68) 
and we have assumed periodic potentials of equal magnitude, but opposite sign, are applied 



on the electron and hole. Analogously to (20), we may transform to the Wannier basis to 



give a electron- hole two-band Hubbard model [TR Pf3] : 



H 



<r n.n' 

4e £ \ u - 



ee [m, n 2 , n 3 , n 4 jeI liCT ,e^ 2CT e Tl 3 (T e Tl4(7 / 
+U hh (ni, n 2 , n 3 , n^ nx<jl h x ^Jx nz(T h n ^i - 2U eh {n x , n 2 , n 3 , n 4 )ej liff 



(69) 
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where 

h 2 



V 2 + (-l) Sih WJr) 



2m, 



Wi(r-n'X) (70) 



U i i(n 1 ,n 2 ,n 3 ,n 4 ) = J d 2 rd 2 r'w*(r - n 1 X)w*(r' - n 2 \)V(\r - r'\) 

x Wj(r' — n 3 \)wi(r — n 4 A), (71) 
U eh (n 1 ,n 2 ,n 3 ,n A ) = J d 2 rd 2 r'w* e (r - (m + (1, l)/2)A)<(r' - n 2 A)V(V|r - r'| 2 + rf 2 ) 
x w ft (r' - n 3 A)w e (r - (n 4 + (1, 1)/2)A), (72) 

with i = e,h and 5^ is a Kronecker delta. The Wannier functions Wi(r) differ for electrons 
and holes due to their different masses. There is a lattice offset of (1, l)/2 since we assume 
a type II potential, i.e. the potential minima locations for electrons and holes differ by half 



a lattice unit. A minimal approximation to (69) is to retain the nearest neighbor terms 



in (70) and on-site terms in (71). Fig. [7] shows the results of our numerical evaluations of 
U = U((n x ,n y ), (n x + l,n y )) = ti((n x ,n y ), (n x ,n v + l)) and Uij = Uij(n,n,n,n). In a similar 
way to Bose- Hubbard parameters of Fig. [5j the application of the periodic potential W e (r) 
acts to increase the electron-electron and hole-hole interaction and decrease the hopping 
amplitude. The electron-hole interaction plateaus off since potential minima of the two 
particle species sit on two spatially separate sublattices. Comparison with Ref. [31] reveals 
that for Uij/(t e + th) ^ 1, and at a density of one exciton per site (half-filling in the 
terminology of Ref. [31]), the excitons will be in a Mott-insulating regime in both the 
electron and hole layers. Thus for a large enough potential Wq the system will lie in such a 
Mott insulating phase. We again assume that the lifetimes of the indirect excitons (which 



can exceed ~ fjs according to Ref. [15] ) should exceed the requirements given in (26) and 



(27) in the respective phases. 

Examining various limits leads us to draw a qualitative phase diagram as shown in Fig. 
8] First consider traveling up the d axis, with Wo = 0. Assuming a periodic potential with 
A = 0.1 yum, one exciton per potential minima corresponds to a density of n exc = 10 10 cm -2 . 
Monte Carlo calculations have predicted that Wigner crystallization should occur at r s pa 37 
[32] 5 corresponding to a density of n wc = 2.3 x 10 8 cm -2 in GaAs. As n cxc is above the Wigner 
crystal melting density n wc , we expect that the bilayer system should be conducting (i.e. 
a metallic phase) for d — > oo. For d — 0, the system is still in a fairly low-density regime 
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(a 2 B /\ 2 <C 1), and thus we expect that the ground state may be described by a BEC (i.e. a 
non-localized metallic state) for sufficiently low temperatures. As d is increased, the Bohr 
radius of the excitons increase, until the exciton wavefunctions start to overlap. Beyond 
this point, the excitons cannot be described as bosons anymore, and the system enters a 
BCS phase [39, 40J. Moving in the direction of increasing Wq for small d, as discussed in 



section IV, we expect a Bose-Hubbard transition into a Mott insulating phase. From the 
considerations of Ref. [41], at unit exciton density we expect the system to be in a Mott 
insulating phase for U ee > U e h, Uhh > U e h, and t4,-/(t e + t h ) ^> 1. We thus expect that a 
transition should occur from the electron-hole plasma phase to a Mott insulating phase for 
large d. Connecting the two boundaries for small and large d leads to the phase diagram 
Fig. |HJ It is plausible to expect that the Bose and Fermi Mott insulating states can be 
smoothly connected, in a similar way to a BEC-BCS crossover [39J SO]- We thus conjecture 
that the first order transition line between the metallic and Mott insulating states can also 
be smoothly connected throughout the phase diagram. The repulsion between the particle 
species generally increases with increasing d, as can be seen in Fig. |4j Thus qualitatively 
the transition should shift to smaller values of Wq for the fermionic limit, as shown in Fig. 

EJ 



VII. SUMMARY AND CONCLUSIONS 



We have considered the effect of applying a periodic potential on interacting exciton- 
polaritons and indirect excitons. Our main result is shown in Fig. [5j where the Bose- 
Hubbard parameters for the on-site interaction U and the tunneling amplitude t was calcu- 
lated. We also derived a guideline (eq. ([58]) ) for the range of parameters necessary to realize 
a phase transition from a BEC phase into a Mott insulating phase. The results suggest 
that very exciton-like polaritons are required to observe the transition. Loosely speaking, 
the reason is that for the typical experimental parameters, the tunneling amplitude t is 
far greater than the interaction energy U . Thus in order to make these parameters on the 
same order, the polariton mass needs to be increased to reduce t. This results in the ne- 
cessity of polaritons with a large exciton component. Alternatively, a very large potential 
amplitude Wq can be applied. The experimental challenge in this case is to maintain U and 
t greater than the experimental temperature and system disorder. Since the energy scale 



18 



of the Hubbard parameters are set by the applied potential period, this favors small A in 
order to increase the energy scale. Although we focused mainly on parameters for GaAs, 
we note our formulas are general enough such that a simple substitution of material and 
geometrical parameters in Fig. [5] should be enough to find the Hubbard parameters for any 
semiconductor system. 

We have also considered the effect of increasing the bilayer separation for indirect excitons, 
where there is a crossover from a Bose-Hubbard model to a double Fermi-Hubbard model. 
The Hubbard parameters for the fermionic limit were derived (Fig. [7]). A Mott transition 
should be present for both limits, thus we argue that there should be a transition for all 
intermediate d. In an analogous way that there is a BEC-BCS crossover for zero potential 
[39| HO], the Mott insulating limit should also crossover from a Bose Mott insulator to a 
double Fermi Mott insulator for large potentials. Our argument is based on connecting 
the various limits of the system, and requires a more rigorous numerical investigation to 
confirm our conjecture. A more detailed investigation of the various phases would require 
an extensive numeric survey of the parameter space, which we leave as future work. 
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Appendix A: Effective polariton interaction 



1. Exciton-exciton interaction 



Eq. (29) may be evaluated by making a change of variables to R = (3 e r e + fihTh and 



p = r e — rh, after which we obtain 

h 2 



U dil (Q,Q',q) 



2E ls + ^ M (Q 2 + Q' 2 ) 



+ 



e iq-Ph{p-p') _|_ e m-Pe{p-p') 



-dq-iqil3 h p+p e p') _ e -dq+iq-{P R p+P hP ')-\ g-A(d)(^/ 'p 2 +d 2 +y l / 'p ,2 +d 2 )/a B _ (^1) 



Eq. (36) may be obtained by performing the p and p' integrals separately and using the 



rotational invariance of q. Figure [2] is obtained by numerically evaluating 
2ir 5 d 

I dil (q,d) = iV£exp(2— X(d)) [l (q, (3 h ,d) 2 + I (q, /3 e ,d) 2 - 2e^I (q, f3 ej d)I (q, f3 h ,d)} , 

qa,B Q>b 

(A2) 



where 



I (q,/3,d)= / drrJ (qa B /3r)exp -X(d)^Jr 2 + (d/a B ) 2 



(A3) 

and Jo{x) is the Bessel function of the first kind. The normalization assuming the electrons 
and holes are confined as delta-functions in the z-direction is 



Nr 



\(dy 



(A4) 



2tt(1 + dX(d)/a B ) 

The electron and hole exchange terms may be obtained by following the derivation given 
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in the Appendix B of Ref. [35]. We obtain (38) and (39), where 

2 



I eKCh (AQ,q,9,/3,d) = (|) N* ex.pl 



(2\(d)d\ 
V a B J 

oo p2tt poo p2tt poo p2tt 



a B J 

POO 

x I dx I d6 x dyi I dQ x I dy 2 I d6 2 xy 1 y 2 
Jo Jo Jo Jo Jo Jo 

x cos{AQa B [(3x cos(6> — 6 X ) + flyi cos(6 — 6>i)] 
+qa B [— xcos9 x — (3yi cos^ + (1 — (3)y 2 cos #2]} 



x exp 



X(d) 



[(y 2 cos 6*2 — y\ cos 9\ — x cos 6 X )' 



+ (y 2 sin 6 2 - y\ sin#i - xsin^) 2 + (d/a B ) 2 } 1/2 

^V* 2 + w« B )' A --' m 



exp 
exp 



2 

\{d) 



exp 



yf + (d/a B ) 2 



yl + (d/a B ) 2 

- v / k \/y{ + x 2 + 2yix cos(#i — a 

1 



\/ 2/1 + x 2 - 2y 2 x cos(6> 2 - X ) \/y% + (d/a B ) 2 \Jy\ + (d/a B ) 2 

(A5) 



2. Saturation interaction 



After substituting the exciton wavefunctions (35) into the expression for the saturation 



interaction (48), we obtain 

Q 



U^(Q,Q',q) 



A 3/2 



-^r ) / d 2 xd 2 r e d 2 r h exp[ — — (\r e -x\ + \r h - x\ + \r e - r h \ 
ira B J J a B 



exp[iQ ■ (f3 e r e + f3 h r h ) + iQ' ■ %x + f3 h r h ) - (Q + Q' - q) ■ (fi e r e + f3 h r h ) - q 

(A6) 

Changing variables to £ = (r e — r h )/(2a B ), rj = ((r e + r h )/2 — x)/a B , and 7 = x/a B , 

UUQ, Q'i v) = J ex P["l^ + v\-\v-t\- 2 \t\] 

exp[ia B £ ■ {p h Q - M + (& - f3 h )q) + ia B r, ■ (q - /3 h Q - f3 e Q% (A7) 
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The dimensionless integral appearing in ( 50 ) is 

iUQ,Q',v) = 51/ d^d V de^r] 

x cos{a B f [&Qcos(0 - - f3 e Q' cos(0' - %) + (& - /3 h )gcos%] 
+ obt? [— cos (^ - - P e Q' cos(6' - 8 V ) + gcos^]} 
x exp 



-\ e + V 2 + 2^ cos(% - 6 V ) - Je + V 2 ~ 2^ cos(% - 0„) - 



where 6 (#') is the angle between qr and Q (Q')- 

Appendix B: Normalization of exciton number states 

By definition, we have 

JV(n) = i(0|(6H6t)»| ). 
n! 

Using the commutation relation (|59]), we obtain 

71—1 

iV( n ) = AT( n - 1) - — ^(0|(6) n - 1 (6 t ) m ^(6 t ) n - 1 - m |0). 

' m=0 

To obtain a simplified expression for the second term, it is useful to define 

[K, 6 f ] = J ] 

where 



f = 2 J d 2 r e d 2 r' e d 2 r h d 2 r' h 4> d (r e - r' h )<f) d (r' e - r' h )<f) d (r' e - r h ) 

xw{P e r e + p h r' h )w{P e r' e + p h r' h )w{P e r' e + (5 h r h )e\r e )h\r h ). 

Retaining only powers linear in the operator K, we obtain 

(0\{tf) n {b) m K{b) n - m \0) = (n-m)n\(0\bKtf\0) 



Substituting this into (B2), we obtain 



N(n) = N(n - 1) - ^(n - l)<0|fe J fiT6 t |0>. 



This definition may be used recursively, to obtain the final result (63) 
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The matrix element above may be calculated according to 

(0\bKtf\0) = I w (W ,d) = 2 J d 2 r e d 2 r h d 2 r' e d 2 r' h Mre - r h )<p d {r e - r' h )MK - r h )M< ~ r' h ) 
xw(fi e r e + (3 h r h )w((3 e r e + (3 h r' h )w((3 e r' e + (3 h r h )w{(3 e r' e + (3 h r' h ). (B7) 

Making the transformation to center of mass coordinates z r = r e — r^, z' r = r e — r' h , 
z r = r' e - r h , z CM = I3 e (r e + r' e )/2 + /3 h (r h + r'J/2, we obtain 

J BP (^o,^)~2 J d 2 Z CM \w(z CM )\ 4 J d 2 Z r d 2 4d 2 4<f) d {z r )M<)M4)M Z r ~ *r + K), (B8) 

where we have used the fact that the relative wavefunction <p{r) extends out to a distance 
of the order of ~ Ob, while the Wannier function extends out to at a distance ~ A, with 
A ^> as- The Wannier function has dimensions of the inverse length (in 2D), hence the 
order of magnitude of the first integral is ~ 1/A 2 . The order of magnitude of the second 
integral is ~ a 2 B , making the whole integral of the order of ~ a 2 B /X 2 . As can be shown by 
direct calculation, integrals involving higher powers in the operator K involve higher powers 



of a 2 B /X 2 . Therefore, the approximation made in (B5) is thus reasonable as long as a 2 B -C A 2 . 
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FIG. 1: (Color online) Schematic device configurations considered in this paper, (a) Exciton- 
polaritons and (b) indirect excitons in a periodic potential. A dc voltage is applied on the indirect 
excitons such that the holes and electrons occupy the top and bottom quantum wells (QW) re- 
spectively. Exciton-polaritons are formed by coupling the excitons to a distributed Bragg reflector 
(DBR). 




FIG. 2: (Color online) The direct exciton-exciton interaction integral I dir (q,d) for three quantum 
well separations d/as = (dashed line, right axis) and d/as = 1,2 (solid lines, left axis). 
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FIG. 3: (Color online) The momentum transfer dependence of the electron-electron exchange 
interaction integral I cxch (AQ = 0, q, 9, f3 e , d) for djas = 0,1,2 (solid lines, left axis) and all 6. 
The photon momentum dependence of the saturation interaction integral J sat (0, 0,q) for djas = 
(dashed line, right axis). 





FIG. 4: (Color online) The direct Coulomb interaction I dir (q = 0,d), the negative of the electron- 
electron exchange interaction — 7 cxch (AQ = 0, q = 0, 9, (3, d), and their sum versus the quantum well 
separation d for all 9 and /3. 
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FIG. 5: (Color online) (a) Dependence of the on-site Coulomb energy JJ coul for three bilayer 
separations d with the periodic potential amplitude Wq. Energies are normalized to units to = 
&m' \a 2 anc ^ ^o° ul = 2e ^x2 ■ 3 ) Dependence of the nearest neighbor hopping t (dashed line, 
right axis) and the on-site saturation energy C/ 8at (solid line, left axis) with the periodic potential 

2 

amplitude Wo. The saturation energy scale is Uq^ = 2hg^J^\u\ 2 'Re{uv*)-^ . For example, for 
polaritons in GaAs with A = 0.5 /im, as = 10 nm, 2hg = 15 meV, and m pol = 0.01M the energy 
scales are t = 0.91meV, U^ oul = 3.58/ieV, and = 0.34//eV. 
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FIG. 6: (Color online) Boundary between Bose and Fermi behavior for the quantum well separation 
d (in units of the Bohr radius ob), the density n and the potential Wq. Well separations less than 
the indicated value indicate bosonic behavior, while Fermi behavior holds for the larger separations. 
"Fermi" behavior means that the fermionic nature of the electron and hole making up the excitons 
become important. 
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FIG. 7: (Color online) Dependence of the nearest neighbor hopping for electrons t e and holes th 
(dashed lines, right axis), as well as the electron-electron on-site Coulomb energy U ee , the hole- 
hole on-site Coulomb energy Uhhi an d the the electron- hole on-site Coulomb energy U e h for A > (i 
(solid lines, left axis). Energies are normalized to units to = gjf^ an d Uq = For example, 
for A = 0.1 /xm, as = 10 nm in GaAs, the energy scales are t = 0.561meV and U F = 3.48meV. 
Parameters for GaAs are assumed here, with m e = 0.067mo, = O.lm-o, and e = 13eo, where mo 
is the free electron mass and eo is the permittivity of free space. 




FIG. 8: (Color online) Qualitative phase diagram of the electron-hole bilayer system at n = 1 unit 
filling density. The labeled phases are BCS phase (BCS), Bose-Einstein condensate (BEC), Bose 
Mott insulator (BMI), and Fermi Mott insulator (FMI). Shading represents the transition from 
boson (light) to fermion (dark) behavior. Solid line denotes a first-order phase transition. 
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